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Abstract 

We suggest that the gauge-invariant hedgehoglike structures in the Wilson loops 
are physically interesting degrees of freedom in the Yang-Mills theory. The trajec¬ 
tories of these “hedgehog loops” are closed curves corresponding to center-valued 
(untraced) Wilson loops and are characterized by the center charge and winding 
number. We show numerically in the SU(2) Yang-Mills theory that the density 
of hedgehog structures in the thermal Wilson-Polyakov line is very sensitive to 
the finite-temperature phase transition. The (additively normalized) hedgehog line 
density behaves like an order parameter: the density is almost independent of the 
temperature in the confinement phase and changes substantially as the system en¬ 
ters the deconfinement phase. In particular, our results suggest that the (static) 
hedgehog lines may be relevant degrees of freedom around the deconfinement tran¬ 
sition and thus affect evolution of the quark-gluon plasma in high-energy heavy-ion 
collisions. 
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1 Introduction 

Topological structures in fields are as important as the fields themselves. 
The well-known realization of this statement is the three-dimensional Georgi- 
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Glashow (GG) model, which has topological defects called’t Hooft-Polyakov 
(HP) monopoles [1]. The presence of these defects determine special features 
of the vacuum structure [2]: external electrically charged particles show linear 
conhnement at large separations while the magnetic fields are screened at large 
distances. Neither the conhnement nor the mass-gap generation are realized 
if the monopoles are absent in the vacuum (like, for example, in a case of the 
noncompact Abelian gauge theory). 

Another textbook example is conventional superconductivity: in the Ginzburg- 
Landau approach the density of the superconducting electrons is described by 
the order-parameter held corresponding to an ehective held of paired elec¬ 
trons (Gooper pairs). The topological structure in the Gooper-pair held is the 
Abrikosov vortex, which drastically inhuences the conduction and thermody¬ 
namic properties of conventional superconductors (for example, subjected to 
an external magnetic held). 

The concept of topological defects in helds was invoked to understand color 
conhnement. This phenomenon, which is one of the most interesting issues 
in QGD, is governed by gluons dynamics described by the pure Yang-Mills 
(YM) theory. The YM theories have zero-dimensional topological objects (in- 
stantons), but these objects cannot explain the color conhnement [3] similarly 
to the GG model [2]. On the other hand, the YM theories lack topologically 
stable monopole-like and vortexlike classical conhgurations. Nevertheless, the 
conhnement was suggested to be driven by condensation or, equivalently, by 
proliferation of the vortexlike [4] and monopole-like [5] nonclassical structures. 
Unfortunately, these structures alone are not supported topologically by the 
SU{N) color symmetry of the YM theory in the continuum space-time. There¬ 
fore, the current understanding of color conhnement is still unsatisfactory. 

On the other hand, YM theory can be entirely formulated [6] in the language 
of colorless composite helds (Wilson loops), which are traces of the loop func¬ 
tionals of the gluon helds. Mathematically, these loops are as good as the gauge 
helds from the standpoint of describing properties (in particular, nonpertur- 
bative) of YM theories [6]. Having the importance of the order-parameter 
(Gooper-pair) held in the Ginzburg-Landau model of superconductivity in 
mind, we compare the properties of this held with the well-known properties 
of the Wilson loop variable in YM theory: 

i. Expectation values: The expectation value of the Gooper-pair held corre¬ 
sponds to the density of superconducting electrons, while the vacuum ex¬ 
pectation value of the Wilson loop in YM theory is related to the interquark 
interaction. Both expectation values are order parameters (nonlocal in the 
case of the Wilson loop): the superconducting phase in ordinary supercon¬ 
ductors is determined by a nonzero value of the Gooper-pair held, while the 
conhning phase in the YM theories is manifested by an area-like asymptotic 
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behavior of the Wilson loop. 

ii. Correlation functions: The correlation function of the Cooper-pair helds 
determines the type and spectrum of the superconducting medium, while 
correlators of the Wilson loops provide information about the glueball spec¬ 
trum of the strongly interacting medium, the intermeson interactions, etc. 
in. Topology: the topological structure in the Cooper-pair held is the celebrated 
Abrikosov vortex. What is an analog of this structure in the gluon helds of 
YM theory in the language of the Wilson loop variable? The answer to this 
question may be important for understanding color conhnement because 
the Wilson loop operator itself is intimately related to this important phe¬ 
nomenon. Moreover, the mentioned structure in the Wilson loop must be a 
gauge-invariant quantity like the conhnement phenomenon. 

General arguments were given in [7] to identify the defect-like trajectories in 
the Yang-Mills theory with the closed paths in space for which the (untraced) 
Wilson loop takes its value in the center of the color group. The main argu¬ 
ment in [7] is that the gluon helds in the vicinity of such “extremal” loops 
are similar to the gauge-invariant hedgehoglike structures, which resemble the 
HP monopoles in the GG model. Below we refer to these structures as to the 
“hedgehog loops”. Developing the idea further, we here investigate properties 
of the hedgehogs loops at a hnite temperature. In particular, the static (ther¬ 
mal) hedgehoglike structures - which are also called below as the “hedgehog 
lines” - should dehnitely arise in an effective theory of Polyakov lines [8], which 
is proposed to describe basic features of the hnite-temperature QCD transi¬ 
tion. Moreover, the thermalization properties of the Polyakov line “conden¬ 
sate” [8] give a qualitatively encouraging insight regarding particle production 
in high-energy heavy-ion collisions [9]. As we show below, the density of the 
hedgehog lines in the deconfinement phase increases rapidly as temperature 
increases. This observation suggests that the hedgehog lines, corresponding to 
the extremal values in the Polyakov line variables, may be important degrees 
of freedom in quark-gluon states created in high-energy heavy-ion collisions. 

The content of this paper is as follows. In Sec. 2, we dehne the hedgehog 
structure in the Wilson loop. Section 3 is devoted to a general discussion of 
hedgehog loop properties. In Sec. 4, we derive an expression for the density 
of the static hedgehogs (the hedgehog lines) as a function of the trace of the 
Polyakov line. We present the results of our numerical simulations in Sec. 5 
and state conclusions in the last section. 


2 Center-valued Wilson loops as hedgehog loops 


The hedgehogs in the Wilson loops are similar to the HP monopoles in the 
GG model [7]. In this model the triplet Higgs held <h = of the static 
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HP monopoles [1] vanishes at the center of the HP monopole x = xq, 


l>(xo) = 0. (1) 

Moreover, in spatial vicinity of a HP monopole the Higgs field has a generic 
hedgehog structure in the in a gauge where the gauge held is regular : 


4(x) = yy“Uo) ■ (I - XaY + 0((x - x„)=) . (2) 

Hedgehog structure (2) inevitably implies that Higgs held (1) vanishes, while 
the converse is false in general case. On the other hand, if the triplet held 
$ vanishes at the point xq, then a general Taylor expansion of the held in 
the vicinity of this point starts from linear terms, and expansion (2) therefore 
corresponds to a generic situation. Note that the matrix Y is assumed to be 
nondegenerate in a general case. 


det F ^ 0 . 


In YM theory,a hedgehog structure can be entirely dehned in terms of Wilson- 
loop variables [7]. We consider an untraced Wilson loop beginning and ending 
at the point xq on the closed loop C, 


Wc{xo) = P exp{i5( dXf, A^{x)'^ . (3) 

c 

This quantity transforms locally as an adjoint operator with respect to the 
SU(2) color symmetry, Hc(xo) —h2^(xo)lTc(xo)ff(xo). To improve the anal¬ 
ogy with the triplet Higgs held <h, we subtract the singlet part from (3) 

tc{xo)=Wc{xo)-l-^TTWc{xo). (4) 

This is a traceless adjoint operator similar to the held <h. 

The next step is to associate the triplet part Pc of Wilson loop (4) with 
the triplet Higgs held <h(a;) in the GG model. To uncover a link between the 
hedgehog structures in the helds of the GG model and in the extremal loops 
of the YM theory, we notice the following [7]: 

i. The Higgs held vanishes in all points x belonging to the HP monopole 
trajectory, Eq. (1). Similarly, one could expect that the Wilson loop eval¬ 
uated on the hedgehog loop C belongs to the center of the gauge group, 
Wc = ±11. In turn, this implies that Pc vanishes on the hedgehog loop C: 
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IVc G Z 2 -v^ Fc — 0. (5) 

We intentionally omit the argument xq in the loop variables in Eq. (5) 
because condition (5) turns out to be independent of the choice of the 
reference point Xq on the trajectory C [7]. In other words, if Wc{xi) G Z 2 
for at least one point Xi G C, then Wc{xq) G Z 2 for all points Xq G C. This 
relation reflects self-consistency of our dehnition of the hedgehog loop. 

In order to illustrate the self-consistency, let us consider two Wilson 
loops lying on the same contour C but open at two different points xq and 
Xi- Both points are arbitrary. Then the corresponding Wilson loops (3) 
are related to each other by the adjoint transformation: 


Wc{xi) = U\xi,Xo) Wc{xo)U{xi,Xo) , 


with 

U{xi,Xo) = P exp^^ig j> <lx^A^,{x)^ G SU{2). 

XQ 

Obviously, if Wc{xo) G ^2 then IFc(a;i) G .^2 as well, and, consequently, 
rc(a;i) = 0. Thus the condition (5) is insensitive to the reference point xq. 

One can alternatively state that the independence of the dehnition of 
the hedgehog on the reference point Xq manifests a conservation of the 
central charge 2 ;, which is dehned in point iv below. 
ii. The triplet helds around the HP monopole have the structure of hedge¬ 
hog (2). To demonstrate that the center-valued Wilson loop variables 
have something to do with the hedgehog structure, we inhnitesimally de¬ 
form the contour C ^ C + 5C in the vicinity of the point Xq- This shift 
also relocates the reference point, Xq ^ x = Xq + Sx (see Fig. 1). Without 



Fig. 1. The Wilson loop C (solid line) and its variation (dotted line). 

loss of generality, we assume that a tangent vector to the loop C at the 
point Xq is aimed in the time direction. A physically distinguishable shift 
Sx should be tangential (spacelike) to the contour C because longitudinal 
(timelike) shifts of the point Xq leave the functional T^ invariant. 

Then the deformed loop rc(a;), in general case, should have a spatial 
hedgehoglike structure in the vicinity of the zero point xq, 
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( 6 ) 


rc(a:) = Y>’?‘(a:o) ■ {x - XoT + 0((x - x„fj . 

This structure is similar to that of Eq. (2) with only one exception: the 
definition of the matrix yc involves the path derivative instead of the 
usual derivative, 


yc{x) = 


6n(x) 


6x^ 


a, f = 1, 2, 3 - 


+X0 


( 7 ) 


which defines a change of the functional T^ under an infinitesimal change 
of the contour C ^ C + 6C. 

It is worth stressing that a loop at which the Wilson loop takes its 
value in the center of the gauge group correspond, in general case, to 
a hedgehog loop. This statement may not be true in a degenerate case, 
dety^^lx) = 0, which is realized, for example, in a trivial vacuum, = 0. 
The degenerate case is nongeneric and thus is statistically suppressed. 

Note that the matrix (7) can also be written in a explicitly Lorentz- 
invariant form, 

yc^'ix) oc , 

where my{x) = Xy/\xy\ is a vector which is tangent to the contour y at 
the point x. The contour y is parameterized by the vector function Xy 
of the variable r and Xy = dxy{T)/dT. 

in. The location xq of the HP monopole is a gauge-invariant quantity. Sim¬ 
ilarly, the trajectory C of the hedgehog loop in the Wilson loop defined 
by condition (5) is gauge invariant because the Wilson loop eigenvalues 
are gauge invariant. 

iv. The SU{2) HP monopoles are classified by the integer-valued winding 
number n G vr 2 (S'^) ~ Z corresponding the number of times the vac¬ 
uum manifold is covered by a single turn around the two-dimensional 
boundary dD^ ~ S'^ of the three-dimensional volume containing the 
HP monopole. Hedgehog structure (2) with oc 5“* corresponds to the 
n=l HP monopole in the GG model. In YM theory, the hedgehog loop 
is classified by the pair of numbers {z,n), where = ±1 is the central 
element of the Wilson loop Wq G Z 2 determined by Eq. (5) and n is 
the winding number associated with hedgehogs structure (6). Below, we 
call 2 : the “center charge” of the hedgehog loop. A hedgehog line, which 
is a particular case of a general YM hedgehog in a Wilson loop, is also 
classified by the two charges. 

V. The gluonic fields around a static hedgehog (or, equivalently, a thermal 
hedgehog line) share similarity with a field of a monopole in the axial 
gauge studied in [10]. However, in general, the hedgehog singularity has 
something to do with electric fields rather than with the magnetic ones. 
Indeed in the case of a static trajectory the matrix y"'^ becomes equal to 
the chromoelectric field, = as one can deduce from Eq. (7) by an 
explicit calculation. 


6 



vi. Because the hedgehog loops are closed, the hedgehog charges must be 
conserved. 


3 Hedgehog loops and color confinement 


As noted in the introduction, a close relation between the hedgehog loops 
and the color confinement problem can be understood intuitively on general 
grounds because the vacuum expectation value of the Wilson loops recognizes 
the color confinement. The hedgehog loops, corresponding to the extremal¬ 
valued Wilson loops, should also be sensitive to the color confinement. In this 
section, we argue in support of this suggestion. 

In conventional superconductivity, Abrikosov vortices are singularities in the 
superconducting condensate (i.e., in the Cooper-pair field). In the core of the 
Abrikosov vortices, the superconductivity is broken, and the normal metal 
state is restored. It is well known that the weaker the Cooper-pair conden¬ 
sate is, the lighter (less tense) the Abrikosov vortices are. As temperature 
increases, the condensate weakens, and nucleation of the vortices due to ther¬ 
mal fluctuations strengthens. Thus, the higher the temperature is, the denser 
the (thermal) vortex ensemble should be. 

It can be expected in the YM theory that the density of hedgehog loops is 
also sensitive to the phase transition. In the Euclidean formulation of the field 
theory, the temperature direction of the system volume is compactified to 
a circle whose circumference is equal to the inverse temperature. The order 
parameter of the phase transition is the vacuum expectation value of the (trace 
of the) Polyakov line. 


l/T 

I/(x) = P expjigf da ;4 A 4 (x, a; 4 )|, (8) 

0 

which is static Wilson loop (3) closed through the boundary of the compact¬ 
ified (temperature) direction (we here define L = ^TrL). Functional (8), also 
called the thermal Wilson line, is a basic variable in an effective theory pro¬ 
posed in [8] to describe the properties of the finite-temperature QCD tran¬ 
sition. In the confinement phase, the expectation value of the Polyakov line 
is zero, {L) = = 0, indicating that the free energy of a single quark 

becomes infinite, Fq oo. In the deconfinement phase, the Polyakov line has 
a nonzero expectation value, {L) ^ 0, and the quarks are no longer confined. 

It can be expected that static hedgehog lines (he., objects corresponding to the 
hedgehog structures in the Polyakov lines) are sensitive to the phase transition 
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because the Polyakov line is the order parameter of this transition. In fact, 
the following analogy between the Cooper-pair condensate and the Polyakov 
line can be established. Both quantities are the order parameters, and the 
singularities in the Cooper pair condensate and the extremal values of the 
Polyakov lines should therefore be physically relevant degrees of freedom in the 
corresponding theories. As we show in the following sections, this suggestion 
is indeed correct. We note that this analogy is only partial because the core of 
the Abrikosov vortex corresponds to the unbroken vacuum while the vacuum 
in the core of the hedgehog loop in the Yang-Mills theory is maximally broken. 


The behavior of the static hedgehog lines can be found qualitatively using sym¬ 
metry arguments. The conhnement-deconhnement phase transition is accom¬ 
panied by breaking the global center symmetry in the deconhnement phase. 
This center symmetry is realized in the case of the SU{2) Yang-Mills theory 
by the transformations L{x) zL{x), where the factor = ±1 belongs to 
the center of the gauge group, G ^ 2 - In the conhnement phase, the center 
symmetry is unbroken, and the effective potential on the Polyakov variable, 
Ves{L) = h 2 L‘^/2 + hiL^/4: +... with 62,4 > 0 , has a minimum at the symmetric 
vacuum (L) = 0, which is far away from the central values L± = ±1. There¬ 
fore, the density of the hedgehogs lines should be small in the conhnement 
phase. 


In the deconhnement phase, the center symmetry is broken by the vacuum, 
and (as the temperature passes the critical value T = T^) the single minimum 
of the ehective potential splits into two ground states, Amin = ±-^— 62/^4 
with 62 < 0, one of which is picked up by the vacuum (A). As the tempera¬ 
ture increases further, the vacuum state approaches one of the central values 
Amin L± = ±1. Therefore, we expect to observe an increase in the density 
of static hedgehog lines characterized by the central charge closest to the 
vacuum state (A). The density of static hedgehog lines with the larger value 
of 1(A) — z\ should decrease in the deconhnement phase. The density of such 
intolerable hedgehog lines for T > T^. should be lower than the density in the 
conhnement phase. It is clear that at the critical temperature, the hedgehog 
line density should split into two branches. 


The qualitative behavior of the spatial Wilson loops is known to be rather in¬ 
sensitive to the deconhnement transition: the area law for the spatial Wilson 
loops is observed in both the conhnement and deconhnement phases. There¬ 
fore, we do not expect a drastic change in the density of spatial hedgehog 
loops as the YM system goes through the phase transition. 



4 Density of static hedgehog lines 


To go beyond the qualitative analysis presented in the preceding section, we 
would like to simulate the YM theory on the lattice numerically. But we imme¬ 
diately notice that the dehnition of the hedgehog loop (see 5 and 6) is rather 
inconvenient from the standpoint of numerical simulations. Indeed, the tra¬ 
jectory of the hedgehog structure in the Wilson loops cannot be determined 
locally as can be done in the case of the Abrikosov vortices in the Ginzburg- 
Landau model of superconductivity or in the case of the HP monopoles in 
the GG model. To check whether a hedgehog loop passes through a particu¬ 
lar point X in the space-time for one particular gauge-field configuration, we 
should analyze the Wilson loops Wc corresponding to all closed loops C that 
pass through the point x 3 C. Obviously, it is difficult to realize this definition 
in numerical simulations of the YM theory on the lattice. 


But at a finite temperature, static hedgehog lines are the most interesting 
because their density should show a noticeable change (in contrast to the 
spatial hedgehog loops) at the phase transition. In a continuum limit, it is 
easy to check whether a static hedgehog line passes through the particular 
point X of the space. For this, we should evaluate (untraced) Polyakov line (8) 
and then check whether it belongs to the center of the group. In the SU{2) 
YM model, the normalized trace of the Polyakov line on a static hedgehog 
line should be equal to ±1. In this section, we derive an analytic expression 
for the (absolute value of the) density of a static hedgehog line as a function 
of the traced Polyakov line p = p[L]. 


Before proceeding further, we point out that in the static limit, hedgehog 
lines are similar to Abelian monopoles in the Polyakov and the Axial Abelian 
gauges of the SU(2) YM theory [11,10,12]. Because the Abelian monopoles 
in the Abelian Polyakov gauge are always static in the continuum limit, they 
are not related to the confinement of static quarks, but they still may be re¬ 
sponsible for the spatial string tension at high temperatures [13]. Moreover, 
such monopoles are related to the topological charge in the YM theory [10,12]. 
The Polyakov-loop variable can also be used to find (static) monopole con¬ 
stituents in physically interesting topologically nontrivial configurations [14]. 
We note that the density of static hedgehog lines (as well as the Polyakov 
line expectation value) is dual in a sense to the Abelian monopole condensate 
because the density of the hedgehog lines is expected to be an order parameter 
of the finite-temperature phase transition while the monopole condensate is a 
disorder parameter [15]. 
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4-1 Density of static hedgehog lines in continuum 


In general, the untraced Polyakov line can be expanded in the quaternion 
basis: 


i/(x, Xi) = -h(x) • 11 + i(f (x, X/fj ■ a ), (9) 

where the singlet part of L is independent of X 4 and the components of the 
vector r describe the gauge-variant triplet part. If a hedgehog line passes 
through the point xq, then P(x) = ±1 and r(xo) = 0. Therefore, the triplet 
part in a neighborhood of the point Xq can be expanded as 


r“(x) = y“*(xo) ■ {x - xqY + 0((x - Xo)^) , (10) 

where det y“^(xo) 7 ^ 0 in general case. We do not consider the rare degenerate 
case dety“^(xo) = 0 here. Therefore the density of the static hedgehog line 
corresponds to the density of the points, where the Polyakov line takes its 
value in the center of the group. We note that Eq. (10) is similar to Eqs. (2) 
and ( 6 ). 

The density of the hedgehog lines is expressed as 


/> = T5(f=) = W-^(|f|), (11) 

where 7 is an unknown factor to be determined below. Because the quantity 
P is dimensionless, the dimension of the factor 7 should be [mass^]. 

Using the identity (5(r) = {r) and Eq. (11), we obtain p[L] = 27r7|f| (5*^^^(f). 

Integrating this identity in the inhnitesimally small volume Uxo with xq 3 
gives the self-consistency condition, which allows us to determine the factor 7 
in Eq. (11): 


1= / dxp(x) = 27r f dx 7 (x)|<h(x)| 5 ®(<h(x)) . (12) 

Changing the integration variables from dx to dP, we obtain the Jacobian 
factor. 


dx 


1 


I det 


ar° I 
dx^ I 


df 


1 -p 

detU“^(xo)| 


(13) 
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Substituting Eq. (13) in Eq. (12), we obtain 


^ _ det _ 1 1 ^^^ (92^2 

^ 27 r|f| ~ 47r-^2(l - L"^) V dx^dx^ ’ 

where we use the identity |r| = vT^-T^ and the relation 


(14) 


d^L\^) 


dx^dx^ 


x=xo 


2y“(xo) E“^'(xo), 


(15) 


which follows from Eq. (10). Noticing that det[2 F^(xo) E(xo)] = 8(dety)^, 
we obtain Eq. (14). Substituting Eq. (14) in Eq. (11), we hnally write the 
expression for the monopole density: 


p[L(x); 


1 1 
16 ^2(1 - L2(x)) 



d2 


dx^dx^ 


L2(x) S fl 

i,j 



( 16 ) 


where E(x) is the Polyakov line evaluated at the point x. The determinant is 
taken over the three spatial indices i,j = 1, 2, 3. 


There are three important properties of formula (16). First, this expression 
depends only on the singlet part (trace) of the Polyakov line and is independent 
on the triplet part. 


Second, Eq. (16) corresponds to the absolute density of hedgehog lines with a 
particular value of the central charge ; 2 . As mentioned, the hedgehog lines are 
described by the winding nnmber n and the center charge ; 2 . Equation (16) 
treats hedgehog lines with ±n winding numbers on eqnal footing while dis¬ 
criminating between different center charges ; 2 . Knowledge of the absolute 
density is enough for our exploratory study reported below. The information 
about the sign of the hedgehog charge is given by the sign of the determinant 
of the matrix K, which can be computed explicitly if the triplet part of the 
(untraced) Polyakov line is calculated. 

Third, Eq. (16) treats the loop variations (7) as ordinary variations because the 
object of our interest is static by dehnition. Equation (16) can be generalized 
to the nonstatic case by properly treating the loop variation discussed in Sec. 2. 


4-2 Density of static hedgehog lines on the lattice 


The density of hedgehog lines in the continnum space-time is given by Eq. (16). 
To obtain the corresponding equation on the lattice, we discretize the double 
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derivatives in Eq. (16) straightforwardly: 


^ /W = /(x - A) - 2/(x) + /(x + /x), (17) 

^2 1 

^^^^/(x) ^ A^A,/(x) = - ^ P^P,/(x + P^-/i + P,-i>), (18) 

where no sum over /i in Eq. (17) is implied. An illustration of the double 
derivatives is provided in Fig. 2. 



h 


Fig. 2. The lattice double derivatives = A^Aj^/(x) at the point x. 

The lattice version of density (16) can be expressed in the form 

P [^(x)] = P+ [i^(x)] + p- [L(x)] , (19) 

where the contributions of the positive and negative center loops to the total 
density can formally be written as 


p± [L(x)] = d± [L(x)] ^(l =F T(x)) . (20) 

The weights of the separate contributions of the positive and negative central 
loops are 


d± [^(x)] 


1 

64 



AiAjL2(x) 


( 21 ) 


Equation (20) is not suitable for use in the lattice simulations. From the practi¬ 
cal standpoint, instead of Eq. (20) it is convenient to evaluate the distributions 


^±(^o) = {d± [T(x)] - T(x))^ . 


( 22 ) 
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Then the expectation values of monopole densities (19) are given by the limits 


(P±) = lim g±{Lo). 

Lo^±l 


(23) 


5 Static hedgehog lines from lattice simulations 


To investigate the static hedgehog lines numerically, we simulated the SU(2) 
gauge theory with the standard Wilson action on the asymmetric {L^ x Lt = 
20^ X 4) lattice with periodic boundary conditions in all directions. We used 
100 independent conhgurations for each value of the gauge coupling, which 
ranges from (3 = 2.2 to j3 = 2.5. The SU(2) gauge theory has the second-order 
phase transition at ~ 2.30 for this lattice geometry [17]. 

Before presenting the results for the center-valued loops, we first discuss some 
well-known properties of the Polyakov lines. The lattice Polyakov line is 


PAU] 


1 

2 


Lt-l 
Tr n 


X4=0 


U4(x, X4) , 


(24) 


where is the lattice gauge held. The distribution of the Polyakov lines. 


D(L) = {S{P,4lU] - L )), (25) 

is symmetric in the conhnement phase and asymmetric in the deconhnement 
phase in a sufficiently large system volume with periodic boundary condi¬ 
tions [16]. These properties rehect unbroken and broken ^2 central symmetry, 
L —L, in the respective conhnement and deconhnement phases. 

Examples of typical distributions obtained for our gluon conhgurations^ are 
shown in Fig. 3. We also show the distribution that is symmetrized under 
the discrete transformation L —L, which may be associated with jumps of 
the system from one ^2 minimum to the other. In the deconhnement phase, 
the symmetric distribution has two maximums, which are known to become 
more pronounced as the volume of the system increases. These maximums 
in distributions correspond to minimums of the ehective potential for the 
Polyakov variables discussed in Sec. 3. 

^ The Polyakov line distributions can be modified by specific boundary conditions 
(for example, twisted) in a finite volume. The properties of Polyakov lines on the 
torus with twisted boundary conditions are discussed in Ref. [18]. 
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(a) (b) 

Fig. 3. Typical distributions (25) of Polyakov lines (24) in the confinement (a) and 
deconfinement (b) phases are shown by open circles. The filled circles correspond 
to L —L symmetrized distributions. 

In the confinement phase, the fluctuations of the Polyakov line are random and 
are therefore basically given by the symmetric Haar measure Zl)Haar(-h)dL = 
(2/7r) Vl — dL. In the deconhnement phase, the distribution is asymmetric 
and, consequently, far from flat. The distributions normalized by the Haar 
measure are shown in Fig. 4. 
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Fig. 4. The same as in Fig. 
D{L) = T)(L)/Vl-T2. 



(b) 

3 but for normalized distributions 


The qualitative behavior of the normalized distributions of the Polyakov line 
are very interesting to us because the Haar factor appears explicitly in the 
dehnition of density of the static hedgehog lines (16). The density of hedgehog 
lines thus explicitly excludes the effect of the Haar measure. The examples in 
Fig. 4 show that the normalized distribution of the Polyakov line is hnite at the 
extremal end L —±1. To obtain an insight into the behavior of the extremal 
ends, we £t the Polyakov distributions by the simple linear function 


D^*(L) = Dz + — z L), as L ^ z = ±1, 


( 26 ) 
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where Dz and Cz are the htting parameters labeled by the center element 2 :. 
Examples of hts (26) are shown in Fig. 5. 




(a) (b) 

Fig. 5. Fits (26) of the normalized distributions D{L) = Zf(L)/\/l — in the 
confinement (a) and deconfinement (b) phases at both sides L ±1 of the distri¬ 
butions. 


We also show the hts of the symmetrized distributions. As can be seen from 
these hgures, the hts work very well. The extreme values of the distributions 
D-i 7 ^ D^i are clearly asymmetric in the deconhnement phase, D±i 7 ^ T*sym = 
^{D_i + D+i) and are symmetric, D±i = in the conhnement phase. 

The extreme values Dz of the normalized distributions D{L) are shown in 
Fig. 6 (a) as functions of the lattice coupling /?. In the conhnement phase 



B B 


(a) (b) 

Fig. 6. (a) The extreme values Dz of the normalized distributions of the Polyakov 
line, (b) The expectation value of the Polyakov line | (L) | as a function of the coupling 
(3. The critical temperature D corresponds to (3c ~ 2.30, and the deconhnement 
phase is realized at /3 > /3c (equivalent to T > Tc). 

{(3 < /5c), the quantities Dz{P) are barely dependent on temperature (/3). But 
as the temperature increases further, the system passes through the second- 
order phase transition, and the extreme quantities Dz{(3) with 2 ; = ±1 show a 
rapid change at the deconhnement point (3 ~ (3c- It can be clearly seen that the 
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values D_i become enhanced while the quantity D+i becomes suppressed in 
the deconhnement phase. This corresponds the system residing in the z = —1 
vacuum ^ . For comparison, the expectation value of the Polyakov line as a 
function of (3 is shown in Fig. 6 (b). 

Distributions (22) corresponding to the contributions q±{L) of the negative 
and positive hedgehog lines (20) to total hedgehog line density (19) are shown 
in Fig. 7 in the conhnement and deconhnement phases. They indicate that the 




Fig. 7. Distributions (22) as functions of the Polyakov line L in the confinement (a) 
and deconhnement (b) phases. 

averaged quantities q±{L) are not symmetric functions of L both in the con¬ 
hnement and in the deconhnement phases: q±{L) ^ Q±{—L). This observation 
stresses the importance of d± factors ( 20 ) in the dehnition of the hedgehog 
line density. It can also be seen that the density distributions in the con¬ 
hnement phase are cross-symmetric in the sense that p±{L) = Qz^{—L). This 
cross-symmetry, which is a direct consequence of the Z 2 center symmetry, is 
violated in the deconhnement phase. 

Finally, the density of static hedgehog lines is given by the extreme values 
L ±1, Eq. (23), of distributions (22). To obtain these densities, we used a 
linear-ht function similar to Eq. (26). The densities p_ and p+ of the hedgehog 
lines with the respective center charges z = —\ and z = +1 are shown in Fig. 8 
as a function of the gauge coupling p. A rapid increase in the density of ; 2 = — 1 
hedgehog lines, which correspond to the value of the center group closest the 
broken vacuum state, can be clearly observed. The hedgehog lines with the 
center charge z = +1 are far from the broken vacuum and are consequently 
suppressed in the inhnitely high temperature limit. The symmetrized density, 
Po = |(P- + P+), shows a less pronounced, but still noticeable, increase as 

^ For illustrative purposes, we always show the cases T)_i(/3) > D_|_i(/3) in the 
deconhnement phase. In a hnite volume, the system jumps from the z = -|-1 vacuum 
to the z = —1 vacuum and back. As the volume of the system increases these jumps 
are suppressed. 
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Fig. 8. The densities pz of static hedgehogs corresponding to the positive and neg¬ 
ative central charges z = +1 and z = — 1 as a function of j3. The vacuum resides at 
the negative center value z = —1. The dotted horizontal line marks the value pconf 
of the densities in the confinement phase. 

the temperature increases in the deconfinement phase. In the conhnement 
phase, all densities take the same value pconf ~ const and are independent 
of the temperature (/3) within error bars. Our results thus indicate that the 
additively normalized density 


Spz(T) = pz(T) - p(T = 0), z = ±l, (27) 

(as well as the corresponding symmetrized quantity 6 po) is an order parameter 
of the deconhnement phase transition. 

Note that in the lattice nnits the density p±. of the hedgehog lines - shown in 
Figure 8 - is more than by an order of magnitude smaller than the extreme 
valnes Dz of the normalized distributions of the Polyakov line, Figure 6(a). 
The difference comes form the d± factors (21), which enter the dehnition of 
the density p± given in Eqs. (22) and (23). 

We have also verihed that the dependence of onr results on the bin size of 
the distribntions (shown in Fig. 7) is very weak. For example, at /? = 2.33 
corresponding to the deconhnement phase, for one of the densities, the valnes 
p+ = 0.0427(17), 0.0455(29), and 0.0468(41), correspond to the respective bin 
numbers Nl = 100, 200, and 300. A similar weak dependence on the nnmber 
of bins is observed for p_. The densities at other values of f3 show the same 
behavior. 

It is important to note the following. One could expect from the very beginning 
that the static hedgehog lines - being dehned by the center elements of the 
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gauge group - should be dense in the high-temperature phase since at very- 
high temperatures the Polyakov line approaches a center element. However, 
this argument is not valid because, for example, it also leads to the erroneous 
conclusion that the hedgehog lines should be dense in the trivial vacnnm, 
= 0. In fact, arguments of Ref. [7] indicate that the trivial vacnnm does 
not possess the hedgehog lines at all. Technically, vanishing of the density of 
the hedgehog lines in the trivial vacnnm is guaranteed by the determinant in 
Eq. (16). Thus, it is not clear ad hoc shonld the density of the hedgehog lines 
be high in the deconhnement phase or not. Conseqnently, the enhancement of 
the density of the hedgehog lines in the high temperatnre phase. Figure 8, is 
a non-trivial dynamical fact. 


6 Conclusions 


We have given analytic and nnmerical argnments supporting the snggestion [7] 
that the gauge-invariant hedgehoglike structnres in Wilson loops (’’the hedge¬ 
hog loops”) are physically interesting degrees of freedom in Yang-Mills theory. 
The hedgehog loops can be dehned self-consistently in terms of the Wilson loop 
variables. The hedgehog loops are associated with the (untraced) Wilson loops 
that take values in the center of the gange gronp, 2 ; G Z^r. In the SU(2) Yang- 
Mills theory, the hedgehog loops are characterized by two quantized numbers: 
the center charge = ±1 G Z 2 and the winding nnmber n G Z. Because the 
hedgehog loops are closed by construction, it seems natnral to snggest that 
the corresponding charges and n are conserved. In fact, the conservation of 
the center charge was proved in [7] and was discussed here. The conservation 
of the winding nnmber n along the hedgehog loop and a stability analysis of 
the hedgehog conhgnrations in particular models may require an additional 
analysis. 

From a general standpoint, it is clear that the hedgehog loops should be re¬ 
lated to the conhnement properties of the system because they are formu¬ 
lated in terms of the Wilson loops, which in turn are intimately related to 
color conhnement in the Yang-Mills theory. The gauge invariance of hedgehog 
loops provides an additional argument in favor of their importance for gauge- 
invariant phenomena, like the color conhnement. Unfortunately, the nonlocal 
formulation of the hedgehogs makes it difficult to locate these objects directly 
in separate conhgnrations of the gluon helds. On the other hand, this difficulty 
can be circumvented by noticing that the statistical properties of hedgehogs 
can be stndied using particular distributions corresponding to these objects 
with trajectories of predehned shapes. These distributions should provide in¬ 
formation about densities, correlation properties, etc., of these objects. 

In this paper, we have demonstrated the success of the proposed numerical 
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approach by studying the properties of static hedgehog lines in the finite- 
temperature SU(2) gauge theory. By definition, static (or, ’’thermal”) hedge¬ 
hog line correspond to extremal (he., center) values of the Polyakov line 
variables, which should naturally arise in an effective model of the Polyakov 
lines [8]. This model is used to describe features of the hnite-temperature QCD 
transition expected to be realized in the high-energy heavy-ion collisions [8,9]. 
We showed numerically that the density of thermal hedgehogs changes rapidly 
as the temperature increases in the deconhnement phase. We observed a sub¬ 
stantial increase (decrease) of the density of the thermal hedgehog lines with 
the central charge closest to (furthest from) the value of the Polyakov line 
in the centrally broken vacuum of the deconhnement phase. We thus found 
evidence that the hedgehog line density (additively normalized to zero at zero 
temperature) is an order parameter of the deconhnement phase transition. 
We conclude that hedgehogs in Wilson loops (in particular, in Polyakov lines) 
may be relevant degrees of freedom of the Yang-Mills vacuum. 
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